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We evaluate the joint distributions of electron waiting times in coherent conductors described by 
scattering theory. Successive electron waiting times in a single-channel conductor are found to be 
correlated due to the fermionic statistics encoded in the many-body state. Our formalism allows us 
also to investigate the waiting times between charge transfer events in different outgoing channels. 
As an application we consider a quantum point contact in a chiral setup with one or both input 
channels biased by either a static or a time-dependent periodic voltage described by Floquet theory. 
The theoretical framework developed here can be applied to a variety of scattering problems and 
can in a straightforward manner be extended to joint distributions of several electron waiting times. 
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I. INTRODUCTION 

Investigations of current fluctuations are useful to 
understand the quantum transport in small electronic 
conductors. 1 In the mesoscopic regime, the transport is 
coherent and the measurement of shot noise provides in¬ 
formation about the effective charge and the quantum 
statistics of the involved quasiparticles. 2-5 A more gen¬ 
eral characterization of the transport is given by the 
full counting statistics (FCS) of transferred charge. 6-9 At 
long times, not only the noise, but all zero-frequency cur¬ 
rent correlators can be obtained from the FCS. Besides 
characterizing the elementary charge transfer events, 10-15 
it has been realized that FCS is also intimately linked to 
entanglement in quantum many-body systems 16-20 and 
to fluctuation theorems at the nano-scale. 21-23 

A complementary picture of the charge transport is 
provided by the electronic waiting time distribution 
(WTD). The WTD is the probability density for a wait¬ 
ing time of duration r to occur between two succes¬ 
sive charges transmitted through a conductor. In recent 
years, WTDs have been considered for various electronic 
systems, including transport through quantum dots gov¬ 
erned by either Markovian 24-28 or non-Markovian 29 mas¬ 
ter equations, and coherent conductors described by scat¬ 
tering theory 30-32 or tight-binding models. 33 For period¬ 
ically driven single-electron sources, 34-36 the WTD pro¬ 
vides a useful characterization of the regularity of the 
emitter and it encodes information about the shape of 
the emitted wave functions which is not readily accessi¬ 
ble in the FCS. 31,3 ' Experimentally, progress has been 
made towards the detection of individual electrons emit¬ 
ted above the Fermi sea, 38,39 possibly paving the way for 
measurements of quantum transport on short time scales. 

The interest in WTDs has so far been focused on the 
distribution of individual electron waiting times. Such 
distributions, however, do not address the question of 
correlations between subsequent waiting times. One may 
ask if the observation of one waiting time will affect the 
following waiting time. Some results indicate that sub¬ 
sequent electron waiting times indeed are correlated. 30,31 
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FIG. 1. (Color online) Joint WTD for a fully open conduction 
channel. The two side panels show the joint WTD along the 
cuts indicated with colored dashed lines in the main panel. 
The joint WTD is well captured by the generalized Wigner- 
Dyson distribution (dashed lines) in Eq. (3). The dotted lines 
in the side panels show the waiting time distributions if there 
were no correlations between subsequent waiting times. Exact 
results have been obtained using Eq. (51). 


However, to fully answer this question, a theory of joint 
WTDs is needed. 

Figure 1 shows an example of a joint WTD. Using a 
method that we develop in this paper, we have calculated 
the joint distribution of waiting times W(ti,T 2 ) between 
electrons transmitted through a fully open conduction 
channel. In Refs. 30 and 32 it was shown that the distri¬ 
bution of individual electron waiting times in this case is 
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well-approximated by the Wigner-Dyson distribution 

Wwd(t) = ^Je- 4 ^ 2 , (1) 
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determined by the applied voltage V. If subsequent wait¬ 
ing times are uncorrelated, the joint WTD in Fig. 1 
should factorize as W(ti,t 2 ) = VVwd(ti)VVwd(t2). 
Such a factorization is often referred to as a renewal 
property. 40 However, as we find, the joint WTD in Fig. 1 
cannot be written in this simple form. This demonstrates 
that subsequent electron waiting times are correlated. In 
fact, based on the analogy between WTDs and level spac¬ 
ing statistics exploited in Refs. 30 and 32, the joint WTD 
is expected to take the form 41 


WWd( t 1j T 2) 
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with b = 729/(1287t). Figure 1 shows that our results 
are well-approximated by this generalized Wigner-Dyson 
distribution, which cannot be factorized into products of 
the WTD in Eq. (1). 

The purpose of this work is to present our method for 
calculating joint WTDs. The method uses a relation be¬ 
tween WTDs and the idle time probability (ITP). 29-31 
The ITP is the probability of detecting no electrons in a 
time interval of length r. For stationary processes, the 
WTD can be expressed as the second derivative of the 
ITP with respect to r. As a second-quantized formu¬ 
lation shows, each derivative corresponds to the detec¬ 
tion of an electron. Building on this principle, we obtain 
WTDs for detections in different channels and joint dis¬ 
tributions of successive electron waiting times. 

To illustrate our formalism, we consider a chiral setup 
where two incoming conduction channels are partitioned 
on a quantum point contact (QPC); see Fig. 2. One or 
both inputs may be voltage biased, either by a static 
voltage or with a series of Lorentzian-shaped voltage 
pulses 1,42,43 as recently realized experimentally. 35,36 We 
calculate the joint WTD and discuss correlations between 
electron waiting times together with WTDs for detec¬ 
tions in different channels. Our findings generalize ear¬ 
lier results for single 30,31 and multi-channel 32 conductors 
and provide insights into the correlations between subse¬ 
quent electron waiting times in phase-coherent conduc¬ 
tors. We focus here on electronic conductors, but our 
concepts may as well be realized in cold fermionic quan¬ 
tum gases. 44 

The article is organized as follows. In Sec. II we discuss 
the theory of WTDs in electronic conductors and explain 
the concept of a renewal process. We calculate the elec¬ 
tronic WTD for the setup in Fig. 2. In Sec. Ill we in¬ 
troduce a generalized ITP which enables us to formulate 



FIG. 2. Chiral setup with electrons in two incoming channels 
being partitioned on a QPC. A (time-dependent) voltage can 
be applied to the contacts of the incoming channels. The 
QPC has transmission T and reflection R. The transmitted 
and reflected electrons are detected in the outgoing channels 
at different positions x = with vf = 1. 


a theory of joint WTDs. In Sec. IV we discuss our re¬ 
sults for the joint WTDs with either individual electrons 
being partitioned on the QPC or electrons from different 
inputs interfering on the QPC in an electronic Hong-Ou- 
Mandel experiment. Finally, in Sec. V we present our 
conclusions. Technical details of our calculations are de¬ 
ferred to two appendixes at the end. 


II. ELECTRON WAITING TIMES 

In this section, we discuss the distribution of waiting 
times between successive electrons emitted into a single 
conduction channel. The WTD is the conditional prob¬ 
ability density of detecting a particle at a time t e given 
that the last detection occurred at the earlier time t s . In 
general, the WTD is a two-time quantity depending on 
both t s and t e . We call a system stationary if the WTD 
depends only on the difference of these times r = t e — f s , 
and we write it as W(r). We denote time intervals by 
t’s, and absolute times by t's. 

We first define the ITP. As we will see, the ITP plays 
a somewhat similar role to that of a generating func¬ 
tion: By differentiating the ITP with respect to its time 
arguments, we obtain the distributions of interest, for 
instance the WTD. This principle becomes particularly 
transparent when using a second-quantized formulation 
of the ITP. We start by deriving the first passage time 
distribution from the ITP before evaluating the WTD. 
We illustrate our formalism by calculating the WTD for 
the setup depicted in Fig. 2 using either a static or a 
periodically modulated voltage. Finally, to gain a better 
understanding of correlated waiting times, we discuss the 
concept of a renewal process. The approach developed in 
this section is important as it allows us to generalize our 
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theory to multiple conduction channels and joint WTDs 
later in this work. 


A. Idle time probability 

The ITP II(t s ,t e ) is the probability of observing no 
transmitted electrons in the time interval [t s ,t e ] at a 
point Xq after the scatterer. To evaluate the ITP we 
employ second-quantization. Working close to the Fermi 
level, we may linearize the dispersion relation as 

e k = hv F k. (4) 

All electrons above the Fermi level thus propagate with 
the Fermi velocity v F along the chiral conduction chan¬ 
nels. We may then evaluate the ITP by considering the 
operator that counts the number of transmitted parti¬ 
cles in the spatial interval x £ [v F t s ,v F t e ]. Here x is 
the distance to the point Xq and the x-axis is oriented 
oppositely to the direction of propagation, see Fig. 2. In 
second-quantization this operator reads 31 



where w(x) and b(x) are the creation and annihilation 
operators of electrons at position x. To keep the notation 
simple, we have omitted the explicit time arguments of 
Q. We will later on make use of the relations 

d t *Q = -&(t 3 )b{t 8 ), 

d t .Q = tf(t e )kt e ), 

setting v F = 1 throughout the rest of the paper. 

The ITP can be expressed as the expectation value of 
the normal-ordered exponential of —<2, 7,31,45,46 

H(O e ) = (: e-« :) , (7) 

where the expectation value is taken with respect to the 
many-body state at t = 0. This state is a non-equilibrium 
state consisting of a Slater determinant of single-particle 
scattering states in the reservoirs. 4 ' To obtain a unidi¬ 
rectional process, we limit our analysis to particles in the 
transport window above the Fermi energy. Depending 
on the specific measurement scheme, the presence of the 
Fermi sea may influence the ITP. This is a subject of 
ongoing investigations. In this work, the notation : • • • : 
denotes the normal-ordering of operators with respect 
to the Fermi sea, 48 i. e., all operators that create elec¬ 
trons above or holes below the Fermi level are moved to 
the left side of the operators that annihilate such exci¬ 
tations, with every individual permutation contributing 
a minus sign. We note that any time evolution of the 
many-body state from t = 0 to t = t° can be absorbed 
into a shift of the spatial interval ( v F = 1) from [f s , t e ] to 
[t a +t°, t e +t°]. With this in mind, we always evaluate the 


expectation value at t = 0 and treat time-dependent sys¬ 
tems by shifting the integration interval in the definition 
of Q in Eq. (5). 

In general, the ITP is a two-time quantity depending 
both on t s and t e . However, for stationary processes it 
depends only on the time difference r = t e —t s such that 

n(t s , t e ) = n (t e - t s ) = n(r), (8) 

just as for the WTD. 

B. First passage time distribution 

To appreciate the usefulness of the ITP, we first show 
how one can obtain the first passage time distribution 
F(t s ,t e ) from the ITP. The first passage time distribution 
is the probability density for the first detection of a parti¬ 
cle to occur at the time t e , given that observations were 
started at the earlier time t s . Such a distribution has 
recently been evaluated for electron transport through a 
quantum dot, with t s being the time at which the exter¬ 
nal electronic reservoirs are connected to the quantum 
dot. 49 To obtain the first passage time distribution from 
the ITP, we notice that the ITP may be expressed as 

U(t s ,t e ) = 1 -J F{t a ,t)dt. (9) 

Here, the integral equals the probability that at least 
one electron is observed in the interval [t s , t e ]. Now, by 
differentiating the ITP with respect to f e , we obtain 

T(t s ,t e ) = -d t .n(t s ,t e )- ( 10 ) 

Furthermore, using Eqs. (6) and (7) we find 

T{t a ,t e ) = (b\t e ) : e"3 : b(t e . (11) 

Equation (11) is very similar to Eq. (7), however, the 
expectation value of the normal-ordered exponential is 
taken with respect to the many-body state with a par¬ 
ticle removed at time t e . This removal constitutes the 
detection event in the definition of the first passage time 
distribution. 

A related quantity is the probability density T(t s ,t e ) 
for detecting a particle at t a with no subsequent detec¬ 
tions until the time t e . Following the line of arguments 
above, we readily find 

T(t a ,t e ) = d t sH(t a ,t e ) (12) 

and 

T(t a ,t e ) = (b\t a ) : e"3 : 6(t s )^ . (13) 

We see that the expectation value of the normal-ordered 
exponential is now taken with respect to the many-body 
state with a particle removed at the initial time t s . For 
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stationary problems, the first passage time only depends 
on the time difference r = t e — t s and we have 

H-r) = Hr) = -drll(T), (14) 

as a direct consequence of Eq. (8) . 

The derivations of Eqs. (10,11,12,13) illustrate an im¬ 
portant principle that we will use in the following: By dif¬ 
ferentiating the ITP with respect to its time arguments, 
pairs of operators are pulled down from the exponential 
function corresponding to detection events at the begin¬ 
ning or at the end of the interval [t s ,t e ]. 

C. Waiting time distribution 

We are now in a position to derive the WTD from the 
ITP. We recall that the WTD is the conditional proba¬ 
bility density of detecting a particle at a time t e given 
that the last detection occurred at the earlier time t s . 
The joint probability density of detecting a particle both 
at t s and t e with no detection events in between is equal 
to the WTD multiplied by the probability density of a 
detection at t s . The joint probability density can be ob¬ 
tained by differentiating the ITP with respect to both 
the initial time t s and the final time t e . Moreover, for 
uni-directional charge transport, the probability density 
of a detection at t s is simply the average particle current 
I[t s ) at time t s . We then find 

I(t s )W{t s ,t e ) = ~d t .d t sll{t s ,t e ), (15) 

where the minus sign comes together with the partial 
derivative with respect to t e , c.f. Eq. (10). Using Eq. (7) 
we now arrive at the second quantized expression 

I(t s )W(t s ,t e ) = (b\t e )b\t s ) : e"* 5 : b(t s )b(t e )) ■ (16) 


and for the WTD 

W(r) = f I(t s )W{t s ,t s + r)dt s . (20) 

' Jo 


D. Applications 

To illustrate our formalism, we consider the setup in 
Fig. 2. Two incoming channels are partitioned on a QPC 
into outgoing channels labeled as a = A,B. The elec¬ 
trons are non-interacting and the temperature is zero. 
The QPC is described by the scattering matrix 



which relates the operators b a (E ) for electrons in the 
outgoing channels to the operators cii(E) for incoming 
electrons originating from contact i = 1,2. The trans¬ 
mission and reflection amplitudes fulfill 

|d| 2 + l r | 2 = T + R = 1 (22) 

as a consequence of current conservation. The incoming 
channels are independently biased with a voltage V l (t) 
which is either constant or periodic in time. For now, we 
keep contact 2 grounded and we focus on the WTDs in 
one of the outgoing conduction channels. The general¬ 
ization to multiple channels and joint WTDs follows in 
the next section. 

The ITP in the outgoing channel a is given by Eq. (7) 
upon attaching an index a to all operators 

n Q (U,t e ) = (:e- < 3“ :}. (23) 


For stationary processes, the average particle current 
equals the inverse mean waiting time, I(t s ) = l/(r). 
Combining Eqs. (8) and (15), we then have 

W(r) = (r)^n(r), (17) 

as previously found in Refs. 30 and 32. For conductors 
driven with period T, we define a one-time ITP by aver¬ 
aging over all possible starting times t s , 31,37 keeping the 
interval t = t e — t s fixed, 


Using the scattering matrix, we may express the Q a 's in 
terms of the operators for the incoming channels. Specif¬ 
ically, we have 

t e 

Qa = R J a\(x)di(x)dx = RQi, (24) 

t * 

and 


n(r) = ^j\(t s ,t s +T)dt s . 


(18) 


This yields the relevant ITP if the observation starts at 
a random time. Employing Eqs. (14, 17), we obtain one¬ 
time distributions which are independent of the detection 
time of the first electron. For the first passage time dis¬ 
tribution, we find 


Hr) 


1 

f 

1 
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[ T{t s ,t s +r)dt s 

Jo 


rT 


T(t s , t s + r)dt s , 


(19) 


Qb=TQ l (25) 

The particles in the incoming channels are in thermal 
equilibrium with the reservoirs. The scattering matrix 
thus allows us to map the evalution of an expectation 
value in a non-equilibrium state onto the evaluation of 
an equilibrium average. Note that the second contact is 
irrelevant since we are only concerned with the transport 
above the Fermi level. 

The single-channel ITP can be evaluated by means of 
the determinant formula 30 , 


o 


Il a (t s ,t e ) =det(I- Q a ) 


(26) 
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FIG. 3. (Color online) Idle time probability (ITP), first passage time distribution, and waiting time distribution (WTD). 
(a) Results for a static voltage applied to contact 1, Vi(t) = V, and no voltage applied to contact 2, V 2 (t) = 0. The QPC 
is tuned to T = 0.3. We consider the electrons that are transmitted into outgoing channel B. The time is given in units of 
t = h/(eV). (b) Results for a series of Lorentzian-shaped voltage pulses of unit charge applied to contact 1, see Eq. (27). The 
width of the pulses is T = 0.02T. The time is given in units of the period of the pulses T. 


where Q a contains the single-particle matrix elements of 
the operator Q a and I is the identity matrix. General 
expressions for a conductor with n channels are given in 
Appendix A. 

We now compare the WTD for a static voltage and that 
of periodic voltage pulses applied to contact 1. With a 
constant voltage, the ITP is a function of the time differ¬ 
ence only and we can directly employ Eq. (17) to obtain 
the WTD. For the time-dependent case, we consider a 
series of Lorentzian-shaped voltage pulses motivated by 
recent experiments. 35,36 For pulses of unit charge, the 
voltage reads (see also Appendix B) 


OO 

eV(t)= ]T 

j = — oo 


2 hr 

(t - jT) 2 + r 2 ’ 


(27) 


where T parameterizes the width of the pulses. Each 
pulse creates a clean single-electron excitation above the 
Fermi level without accompanying electron-hole pairs. 
Such excitations were first proposed theoretically by Lev¬ 
itov and co-workers'’ 42,43 and later named levitons fol¬ 
lowing their experimental realization. 35,36 Importantly, a 
time-dependent voltage has the effect of adding a phase 
to the single-particle states in the contact. 50 We can treat 
this additional phase as a scattering phase which is picked 
up after the particles leave the contact but before they 
arrive at the scatterer, i. e. 

S Vi {t) = e~ i<k St t o dt ' Vi{t ' ) . (28) 


Here t 0 is the time when the voltage is switched on. In 
this way, the contacts can be treated as equilibrium reser¬ 
voirs at the same chemical potential and all the effects 
due to the time-dependent driving are captured by the 
phase Eq. (28). In the energy basis, the effect of the pe¬ 
riodic voltage is that particles can absorb or emit energy 


quanta of size Ml, where Q = 27r/T is the frequency of 
the driving. The Floquet scattering matrix 51 encodes the 
quantum mechanical amplitudes for these processes and 
relates operators before and after the scattering phase as 

OO 

ai(E)= SnK(E-n)- (29) 

n=— oo 

The operator a[ (E) annihilates an equilibrium particle in 
contact i, while a,i(E ) annihilates a particle incident on 
the scatterer. We have also defined the energies 

E n = E + nHQ. (30) 

The S n ’ s are obtained by Fourier transforming Eq. (28). 
Based on Eq. (21), the a,i(E) operators can be related 
to the b a (E) operators and the ITP can be calculated 
analogously to the static case using Eq. (26) (see also 
Appendix A). Finally, to evaluate the WTD, we employ 
Eqs. (17, 18). 

Figure 3 shows the ITP, the first passage time distribu¬ 
tion, and the WTD for the static voltage as well as for the 
train of voltage pulses. Interestingly, the two ITPs are 
nearly indistinguishable despite the very different volt¬ 
ages applied to contact 1. In both cases, the ITP decays 
monotonously with time from n(0) = 1 at r = 0. Turn¬ 
ing next to the first passage time distribution, some struc¬ 
ture starts to be visible and the two voltage cases can now 
be distinguished. Finally, considering the WTDs, we see 
how they clearly differentiate between the two voltages. 
For the static voltage, the electronic wave packets are 
strongly overlapping and the WTD is rather structure¬ 
less. By contrast, the application of Lorentzian-shaped 
voltage pulses leads to the emission of well-localized elec¬ 
tron wave packets with a small overlap, giving raise to 
clear peaks at multiples of the period of the driving T. 
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For a fully transmitting QPC (T = 1), the WTD essen¬ 
tially consists of a single peak around r ~ T (see also 
Ref. 31). For a QPC with a finite transmission, as in 
Fig. 3, levitons may reflect back on the QPC, giving rise 
to the series of peaks in the WTD. We note that the 
period and the width of the Lorentzian-shaped pulses, 
which determine the overall structure of the WTD, are 
easily tunable in an experiment. 

Independently of the applied voltage, the WTD is sup¬ 
pressed to zero at r = 0. This is a consequence of the 
Pauli principle which forbids two electrons to occupy the 
same state. For a QPC with non unit transmission T, 
the mean waiting time takes the form 


(t)t = 


( t )t =i 

T 


(31) 


where the value at full transmission, (t)t= i, is given by t 
defined in Eq. (2) for the static voltage and by the period 
of the driving T for the levitons. 


E. Renewal theory 

It is instructive to investigate the influence of the QPC 
in more detail. In this connection we also introduce the 
concept of a renewal process. The effect of the QPC can 
be understood by resolving the WTD with respect to the 
number of reflections on the QPC that have occurred. 31 
Specifically, we expand the WTD as 

OO 

W(r) =T^R"W ( ” ) (r), (32) 

n —0 

where W^ir) is the WTD between transmitted elec¬ 
trons given that n reflections have occurred during the 
waiting time r. These n-resolved WTD’s can be related 
to the joint probability distributions Wi n (Ti,... ,r„) for 
n successive waiting times between incoming electrons, 
corresponding to the joint WTDs at full transmission. 
We have for example 


W (0) (r) = Win(r) (33) 


distributed. 411 The joint probability distributions then 
factorize as 


Win(T 1 ,...,r n )w JJW^Ti). (35) 

i—1 

Consequently, in Laplace space, we have 

W^iz) = W m (z) (36) 

and 

W w (z) = [W in (z)] 2 , (37) 

and generally 

W^ n) {z) = [W in (z)] n+1 , (38) 


making use of the fact the n-resolved WTDs under the 
renewal assumption are convolutions in the time domain. 

Laplace transforming Eq. (32), we now obtain a geo¬ 
metric series which can be resummed as 31 


W(z) =T^2R n [W in (z)] n+1 

n—0 

TWin(z) 

l-RW in (z)' 


(39) 


Returning to the time domain using an inverse Laplace 
transform, this expression provides us with a direct test 
of the renewal assumption of uncorrelated waiting times. 
We recall that Wi n ('r) is the WTD at full transmission. 

Figure 4 shows a comparison between exact results for 
the WTD and the approximation based on renewal the¬ 
ory. Both for a static voltage and for levitons generated 
by periodic voltage pulses, the renewal assumption cap¬ 
tures the main qualitative features of the WTDs. How¬ 
ever, upon closer inspection it becomes clear that the 
agreement is not perfect. This is an indication that the 
renewal assumption is not correct and that subsequent 
electron waiting times are correlated. With this in mind, 
we go on to develop a theory of joint WTDs. 


III. GENERALIZED IDLE TIME PROBABILITY 


and 

W (1) (t) = f W in (r',r - r')dr\ (34) 

Jo 

where we have integrated over the time at which the re¬ 
flection occurs. Similar expressions can be formulated 
for the ?r-resolved WTDs of higher orders in terms of the 
joint probability distributions for the incoming electrons. 

In the following section we extend our formalism to 
calculations of joint WTDs. However, at this point 
we can proceed using a renewal assumption: We make 
the assumption that successive waiting times are uncor¬ 
related and thus statistically independent and equally 


We now generalize the concepts from the last section 
to an arbitrary number of channels and to an arbitrary 
number of successive waiting times. To this end, we in¬ 
troduce a generalized ITP from which we obtain the joint 
WTDs as well as other distributions. In the following sec¬ 
tion we illustrate our formalism with specific examples. 

We consider a scatterer which is connected to Nj in¬ 
coming and N a outgoing channels. The generalized ITP 
is the probability that no particles are detected in any 
of the outgoing channels during the channel-dependent 
time intervals [£„,£„]. The generalized ITP reads 

n(£;,£!;...;£^,^ 0 ) = (:e-S-*«“ :) , (40) 
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FIG. 4. (Color online) Comparison with renewal theory. Exact results (solid lines) for two different transmissions of the QPC 
are shown together with the approximation in Eq. (39) (dashed lines) based on a renewal assumption of uncorrelated waiting 
times, (a) Results for a static voltage applied to contact 1, Vi (t) = V, and no voltage applied to contact 2, V 2 (t) = 0. The 
QPC is tuned to T = 0.3 (green lines) or T = 0.7 (blue lines). We consider the electrons that are transmitted into outgoing 
channel B. The time is given in units of (r) = h/(TeV). (b) Results for a series of Lorentzian-shaped voltage pulses of unit 
charge applied to contact 1. The width of the pulses is T = 0.0277 The time is given in units of the period of the pulses T. 


having defined projectors for each channel 



The evaluation of the generalized ITP is discussed in Ap¬ 
pendix A. The introduction of individual starting and 
ending times for each channel allows us to compute a 
variety of WTDs. In each channel, the idle time inter¬ 
val [f* , t e a ] can be modified and detection events can be 
inserted by differentiation with respect to the time argu¬ 
ments. The single-channel ITP can always be recovered 
by letting the length of the intervals in all other channels 
go to zero 

n(C£) = n (42) 

where the curly brackets imply V i. The operators Q a 
in Eq. (41) count the number of particles in the spatial 
intervals For this reason, the generalized ITP is 

closely related to the joint particle number statistics in 
spatial subregions. s2,53 

A. Multiple channels 

For the sake of simplicity, we consider only two out¬ 
going channels labeled as A and B. The generaliza¬ 
tion to more channels is straightforward. With two 
outgoing channels, the ITP has four time arguments 
n *!?>*!?)• We consider two different types of 

WTDs: The two-channel WTD is the distribution of 
waiting times between detections in either of the two 
channels. The cross-channel WTD is the distribution 
of waiting times between a detection in one channel and 


the next detection in the other channel. These WTDs 
generally have two time arguments but become one-time 
quantities for driven systems by averaging over a period 
following Eq. (20). 

We first discuss the two-channel WTD. This is the con¬ 
ditional probability density of detecting an electron at 
time t e in either channel, given that the last prior detec¬ 
tion happened at the earlier time t s in any of the two 
channels. This WTD follows from the generalized ITP 
by differentiation with respect to the time arguments 

I(t s )W AB (t s ,t e ) = -dtsdt.n (t s ,t e -,t s ,t e ), (43) 

where we have set t s = t s A = t s B and t e = t e A = t e B , since 
we do not differentiate between the two channels, and 
I = I a + Ib is the sum of the particle currents in each 
channel. If the two channels are uncorrelated, the ITP 
factorizes as 32 

n uc (t s A , t e A] f B , t%) = n A (t% t e A )n B (t s B , t e B ), (44) 

where IIo, is the ITP in channel a and the superscript 
uc stands for uncorrelated. In this case, the two-channel 
WTD takes on a particularly illuminating form 

i(t s )WT B (t s ,t e ) = i A (t s )w A (t s ,t e )n B (t s ,t e ) 

+ T B {t s ,t e )T A {t s ,t e ) + A o B. 

The first term represents contributions where both de¬ 
tections happen in channel A, while no detections occur 
in channel B. The second term corresponds to contribu¬ 
tions where the first detection happens in channel A at 
time t s and the next detection occurs in channel B at 
time t e . Finally, the term A •<->• B indicates that the roles 
of the two channels can be interchanged. 
















In contrast to the single-channel WTD, the Pauli prin¬ 
ciple does not force the two-channel WTD to vanish at 
short times (detections in the two channels can occur si¬ 
multaneously). Evaluating Eq. (43) at t s = f e , we find 


W AB {t s ,t s ) = 2 


m 


(46) 


showing that the two-channel WTD at t s = t e is (twice) 
the coincidence rate in the two channels divided by the 
total particle current at t s . The factor of two accounts 
for the fact that either of the two detections can be in¬ 
terpreted as being the first one. For two uncorrelated 
channels, the expectation value factorizes and we find 

wr B (t s ,n = 2 IAts l ) ^ B ) [ts) , ( 47 ) 


where I a is the average particle current in channel a. 

Next we discuss the cross-channel WTD. This is the 
conditional probability density of waiting until t e for the 
first detection in channel B to happen after a detection 
has occurred in channel A at the earlier time t s . We note 
that additional detection may occur in channel A after 
the time t s . The cross-channel WTD follows from the 
ITP as 


I A {t s )W A ^ B {t s X) = - d t s A d t Al(t s A ,t s -,t s ,t e ) 


= (tf A (t s )bUn ■■ :b B (t e )b A (t s )). 

(48) 

We set t A = t A = t s after taking the derivatives since 
additional detections in channel A may occur following 
the first detection. For two uncorrelated channels, we 
recover the first passage time distribution in channel B 
since the detection in channel A only defines t s without 
influencing channel B 


W^ B (t s ,f) = J s (t s ,f). (49) 


For t s = t e , we again obtain the coincidence rate in 
the two channels. However, this time without the factor 
of two since the event in channel A by definition is the 
first one 


from the multichannel ITP defined in Eq. (40) by intro¬ 
ducing auxiliary channels. For the sake of brevity, we 
consider the case of just two successive waiting times in 
a single channel. The extension to several channels or 
waiting times is straightforward. 

To find the joint WTD, we consider the two-channel 
ITP n(f* , t e a ; t s p, fg), where /3 denotes an auxiliary chan¬ 
nel. Eventually, we set tp = tp = t m and a = /3 and 
skip the channel index. Specifically, we express the joint 
WTD as 


i(t s )w{t s ,t m A e ) = d t ' a dt % d t% n(cc^)| ts=te=t 


= (tf {t m )b ] (t e ) : e~ Q : b(t e )b(t m )b(t s )^ , 

(5i) 

where the operator Q is given in Eq. (5). Using the aux¬ 
iliary channel a detection event is inserted at time t m 
between the starting t s and the end time t e of the inter¬ 
val [t s , t e ] . In a similar way, additional detection events 
within the interval can be introduced. Equation (51) is 
a central result for the joint distribution of two succes¬ 
sive electron waiting times in a single conduction channel. 
Based on this expression we calculated the joint WTD in 
Fig. 1 and we will make further use of it in the follow¬ 
ing section, where we illustrate our method with specific 
examples. 

For a driven conductor, a two-time WTD is obtained 
by integrating over the period T, cf. Eq. (20), 


kV(r 1 ,T 2 ) 


Y £ +n,t s + r 1 + r 2 )dt s . 


In addition, we recover the standard WTD by integrating 
over the time at which the last detection event occurred 


OO 

W{t s ,t m ) = J dt e W(t s ,t m ,t e )- (52) 


As discussed in Sec. HE, the joint WTD appears in 
the expansion of the WTD for a QPC. Returning to the 
results in Fig. 4, the ITP reads 

n(U,t e ) = e ~ T ^ 1 ) (53) 


W A ^ B (t s ,t s ) 


( UnW )) 

W) 


(50) 


For uncorrelated channels this expression reduces to 
I B (t s ) in accordance with Eq. (49). 

In the following section, we evaluate the two WTDs 
and show how waiting times may be correlated. 


B. Multiple times 

We are now in position to formulate our theory of joint 
WTDs. The joint distribution of n waiting times is the 
conditional probability to find a given sequence of n wait¬ 
ing times. As we will see, the joint WTD can be obtained 


where Q\ in Eq. (24) acts on the incoming channel 1. We 
can formally expand this expression in either T or I?. 30,31 
To zeroth order in T, the ITP is unity as no electrons are 
transmitted through the QPC. The n’th order term in 
the expansion yields the reduction of the ITP due to the 
probability that n particles were transmitted through the 
QPC. To zeroth order in R, the ITP equals the ITP for 
a fully transmitting QPC. The n’th order term in the 
expansion in R equals the increase in the ITP due to the 
probability that n particles were reflected. 

Here we perform an expansion in R , and by differenti¬ 
ating the ITP with respect to t s and t e we obtain 

OO 

W(t s ,t e ) =T^2 R n yV (n) (t s ,t e ) (54) 

n =0 
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FIG. 5. (Color online) Joint waiting time distributions for a static voltage, (a) Results for the voltage Vi(t) = V applied to 
contact 1 and no voltage applied to contact 2, Vi (t) = 0. The QPC is fully transmitting T = 1 and we consider the electrons 
that are transmitted into outgoing channel B. The time is given in units of (r) = h/(TeV). (b) Difference AW(ti,T 2 ) = 
W(ti, T 2 ) — W(ti)W(t 2 ) between exact results and results for uncorrelated waiting times. Positive correlations are indicated 
with red, while areas of negative correlations are blue. (c,d) Similar results for a half-transmitting QPC, T = 1/2. 


where 


t e 

W (n \t s ,t e ) = /dt 1 • • • di n W in (tW ■ • ,iV e ) 

n! J 

(55) 

is the WTD given that n reflections on the QPC have 
occurred and W m (t s , t 1 , ■ ■ ■ ,t n ,t e ) is the joint WTD for 
n + 2 detection events in the incoming channel 1. Each 
term in the expansion is the probability density for n par¬ 
ticles to be reflected (corresponding to the prefactor R n ) 
followed by a transmission (corresponding to the pref¬ 
actor T). We integrate over all possible times that the 
reflections can occur and the factor 1/ro! corrects for mul¬ 
tiple counting of reflections. We see that joint WTDs 
occur already in the expansion of the WTD for a QPC. 
By averaging over t s and making a renewal assumption, 
we can recover the results from Sec. HE. 


IV. CORRELATED WAITING TIMES 

We now illustrate our formalism for joint WTDs using 
the setup in Fig. 2. First we consider the partitioning of 
electrons emitted from one source. In the second example 
we discuss an electronic analog of the Hong-Ou-Mandel 
interferometer, where electrons from different input chan¬ 
nels interfere on the QPC. 54-59 


A. Single-source partitioning 

We consider the setup in Fig. 2 with contact 2 
grounded and contact 1 biased with a constant voltage 
or a train of Lorentzian-shaped voltage pulses. Figures 3 
and 4 show single-channel WTDs for this setup. We now 
go on to calculate two-channel and joint WTDs. 

The two-channel WTD is the distribution of waiting 
times between successive electrons irrespective of the 
channel in which they are detected. As the QPC merely 
distributes the incoming electrons into the two outgoing 
channels, we expect simply to recover the WTD of the 
incoming electrons from channel 1. Indeed, using our 
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FIG. 6. (Color online) Joint waiting time distribution for levitons. (a) Results for a series of Lorentzian-shaped voltage pulses 
of unit charge applied to contact 1. The width of the pulses is T = 0.057” and the QPC is tuned to T = 1/2. The time is given 
in units of the period of the pulses T. (b) Difference A>V(ti,T 2 ) = W(ti,T 2 ) — W(ti)W(t 2 ) between exact results and results 
for uncorrelated waiting times. Positive correlations are indicated with red, while areas of negative correlations are blue. 


formalism we find from Eq. (43) that 


within our formalism we readily find, cf. Eq. (48), 


I(t s )W AB (t s ,t e ) = -d t sd t e (: 

= -Me (: e~& 


(56) 


where the integrations in the Q operators run from t s 
to t e and I(t s ) is the average current in channel 1. In 
addition, we have used Eqs. (24) and (25) to relate Q a 
to Q i combined with Eq. (22). For a constant voltage, 
the WTD was evaluated in Ref. 30 and found to be well- 
approximated by a Wigner-Dyson distribution reflecting 
Fermi correlations between the incoming electrons. For a 
train of Lorentzian-shaped voltage pulses, the WTD was 
evaluated in Ref. 31 and found to be peaked around the 
period of the driving with small satellite peaks due to 
the finite overlap of the voltage pulses. As we see here, 
the QPC has no effect on the WTD if we consider de¬ 
tections irrespective of the channel where they happen. 
This argument only applies if the QPC partitions elec¬ 
trons emitted by a single source. In that case, detections 
in the outgoing channels are not statistically indepen¬ 
dent, and the suppression of the WTD at t = 0 remains. 

Next, we consider the cross-channel WTD. This is the 
distribution of waiting times between a detection in one 
outgoing channel and the next detection in the other 
channel. Since the QPC just randomly partitions the 
incoming electrons into the outgoing channels, we ex¬ 
pect the cross-channel WTD to equal the single-channel 
WTD. The cross-channel WTD is conditioned on the first 
detection happening in channel A , whereas the single¬ 
channel WTD is conditioned on the first detection hap¬ 
pening in channel B. In either case, the reflection or 
transmission of the first particle does not influence the 
particles that traverse the QPC at a later time. Indeed, 


w A ^ B {t s ,t e ) = w B {t s ,t e ) 

= (a\{t s )a\(t e ) : e'^ 1 : ai(f e )ai(C)) , (57) 

having used I A {t s )/R = I B (t s )/T = I(t s ). 

For the partitioning of electrons emitted from a sin¬ 
gle source, we see that the two-channel WTD and the 
cross-channel WTD do not provide additional informa¬ 
tion compared with the WTD itself. This is a direct con¬ 
sequence of the operator proportionality Q A oc Q B and 
it does not hold when both contacts emit electrons above 
the Fermi sea as we will see in the following example. 

We now consider just one of the outgoing channels and 
calculate the joint probability distribution of finding two 
successive waiting times t± and r 2 using Eq. (51). For a 
constant voltage V, the joint WTD at full transmission 
is shown in Fig. 1. As mentioned in the introduction, 
the joint WTD is well-approximated by the generalized 
Wigner-Dyson distribution in Eq. (3). We find that the 
joint WTD is symmetric with respect to an exchange of 
the waiting times, t ± f)T2. This symmetry implies that 
the WTD does not change if we invert the spatial argu¬ 
ments of all the b operators in Eq. (51). The symmetry 
in the WTD is thus a consequence of the spatial inversion 
symmetry of the many-body wavefunction. 

In Fig. 5, we show joint WTDs both for full trans¬ 
mission and at half transmission. To highlight possible 
correlations, we also show the difference 

AW(n, r 2 ) = w(n, r 2 ) - W(ti)W(t 2 ), (58) 

between the joint WTD and a factorized WTD cor¬ 
responding to uncorrelated waiting times. 40 The figure 
clearly demonstrates that electron waiting times are cor¬ 
related. The probability of observing a waiting time 
which is shorter (longer) than the mean waiting time 
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(t) is reduced if the previous waiting time was already 
shorter (longer) than the mean waiting time. On the 
other hand, a short waiting time will likely be followed 
by a long waiting time and vice versa. This conclusion 
holds for both values of the transmission. 

A similar analysis can be carried out for levitons emit¬ 
ted by a train of Lorentzian-shaped voltage pulses. In 
Fig. 6 we show W(r 1 , r 2 ) and AW’(r 1 , r 2 ) for levitons inci¬ 
dent on a QPC tuned to half transmission. Again, we ob¬ 
serve a symmetry under the interchange of waiting times 
due to the spatial inversion symmetry of the fermionic 
many-body state. The joint WTD displays a lattice- 
like structure in contrast to the results for the dc-biased 
source in Fig. 5. This reflects the regular on-demand na¬ 
ture of the leviton source. Also here, the probability to 
observe two short or two long waiting times after each 
other is suppressed, while a short (long) waiting time is 
likely followed by a long (short) waiting time. 

Due to the external driving, the sum of two subse¬ 
quent waiting times is likely to equal a multiple of the 
period. If a particular waiting time is shorter than the 
average waiting time, this will not affect the absolute 
emission time of the next electron. Consequently, the 
next waiting time will most likely be longer. Thus, a 
strong regularity in the absolute electron emission times 
leads to strong dependences of the waiting times on each 
other. For the case that the voltage pulses overlap more 
and more until the limit of a constant voltage is reached, 
the regularity in the emission times is provided by the 
Pauli principle. Moreover, this effect is independent of 
the QPC transmission. In the low transmission regime 
(not shown), however, electron emission becomes increas¬ 
ingly rare and transport resembles a Poisson process. 30 
Correlations between waiting times will then eventually 
be negligible. 


B. Hong-Ou-Mandel interferometry 

As a second application, we consider the electronic 
analog of the Hong-Ou-Mandel experiment developed in 
quantum optics. 54-59 In this case, a driven on-demand 
source in each incoming channel emits single electrons 
onto the QPC in Fig. 2 with a possible time delay At 
between the emissions from the two sources. Such exper¬ 
iments have recently been realized with two mesoscopic 
capacitors 59 and two leviton sources. 35 Here we treat the 
periodic emission of levitons onto the QPC. 

Figure 7a shows the distribution of waiting times be¬ 
tween detection events in either of the two outgoing chan¬ 
nels. With zero time delay, incoming levitons antibunch 
on the QPC and there is a large probability of detect¬ 
ing the two outgoing levitons nearly simultaneously. The 
peak in the WTD around the period T corresponds to 
the waiting time between the last leviton in one cycle 
and the first one in the next cycle. Since the probabil¬ 
ity of measuring the waiting time between two levitons 
within one cycle is equal to the probability of measuring 


the waiting time between levitons in different cycles, the 
areas underneath the two peaks are the same. 

In general, for a given time delay At, we find peaks 
in the WTD at At and at T — At, corresponding to the 
waiting times between levitons within the same cycle and 
the waiting between levitons in different cycles. In the 
special case At = T/2, the two peaks merge into one, 
except for a small remaining feature at t = 77 This fea¬ 
ture is related to the satellite peak seen in the WTD for 
levitons in a single channel with perfect transmission. 31 
The WTD decays strongly for waiting times beyond the 
period, since all levitons will be detected independently 
of the QPC. This reflects the high reliability of the leviton 
sources which emit exactly one electron per cycle. This 
example demonstrates how the two-channel WTD pro¬ 
vides information about the synchronization of the two 
sources. Importantly, this detailed characterization does 
not depend on the transmission of the QPC. 

In Fig. 7b we turn to the distribution of waiting times 
between detections in different channels. This distribu¬ 
tion is rather similar to the WTD for a single channel, 
shown in the left inset for comparison. Unlike the single¬ 
channel WTD, the cross-channel WTD is not suppressed 
to zero at t = 0, since detections can occur simultane¬ 
ously in the two channels. For finite delay times between 
the two sources, the particles can go into the same out¬ 
going channel, and both WTDs show a peaked structure 
even for waiting times larger than the period. 

Due to the Pauli principle, two levitons arriving si¬ 
multaneously at the QPC will go into different outgo¬ 
ing channels with unit probability, independently of the 
transmission. Thus, for zero time delay all results are in¬ 
dependent of the QPC transmission T and reduce to the 
results for full transmission. As the time delay between 
the two sources is increased, the overlap between levitons 
arriving at the QPC decreases, and simultaneous detec¬ 
tions in the outgoing channels become increasingly rare. 
For At = T/4, the cross-channel WTD essentially shows 
four peaks within the first period. The peak at t = 0 
is a relict of the fermionic antibunching due to the finite 
overlap of the pulses. The peaks at At and at T — At 
are caused by two successive particles entering opposite 
arms of the interferometer, while the peak at T corre¬ 
sponds to two successive particles entering the same arm 
and the next particle entering the opposite arm. For the 
maximal detuning At = T/2, the WTD resembles that 
of just a single channel. 

Due to the Pauli principle the cross-channel WTD of 
a single source is suppressed for zero waiting times. By 
constrast, with two sources we observe a local maximum 
due to the antibunching of electrons that arrive simulta¬ 
neously at the QPC. The right inset of Fig. 7b shows the 
cross-channel WTD at t = 0 as a function of the time de¬ 
lay, with a maximum for zero time delay in analogy with 
the Hong-Ou-Mandel peak found in the zero-frequency 
current cross-correlations. 58,59 

We see that the two-channel and cross-channel WTDs 
contain different information about the emission pro- 
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FIG. 7. (Color online) Distributions of waiting times for a Hong-Ou-Mandel experiment with levitons. Two sources emit 
levitons of pulse width V = 0.05 T toward a QPC with transmission T = 1/2. The tunable time delay between the sources is 
denoted as At. (a) WTD for waiting times between events occurring in either of the two outgoing channels for different values 
of the time delay, (b) WTD for waiting times between detection events in different channels. The left inset shows the WTD 
for just a single channel. The right inset shows the WTD at r = 0 as a function of the time delay At between the sources. 


cesses. The two-channel WTD mostly contains informa¬ 
tion about the sources alone, in particular their synchro¬ 
nization. The cross-channel WTD contains additional 
information about the QPC and shows more prominent 
signatures of fermionic antibunching. 


V. CONCLUSIONS 

We have developed a general framework for calculat¬ 
ing joint WTDs in electronic multi-channel conductors. 
The central building block of our formalism is the gener¬ 
alized idle time probability, i. e. the joint probability for 
no detections to occur in the outgoing channels. By cal¬ 
culating the joint WTD for a single conduction channel, 
we have explicitly demonstrated that the electron wait¬ 
ing times in coherent conductors are correlated due to 
the fermionic statistics encoded in the many-body state. 

Drawing on the analogy between random matrices 
and free fermions, we have shown that the joint WTD 
for a fully transmitting conduction channel is well- 
approximated by a generalized Wigner-Dyson distribu¬ 
tion. In contrast to a renewal process with uncorrelated 
waiting times, we find that the probability of observ¬ 
ing a long (short) waiting time following a short (long) 
waiting time is increased, while finding two long or two 
short waiting times in succession is less likely. This holds 
both for electrons coming from a dc-biased contact and 
for levitons emitted on top of the Fermi sea by applying 
Lorentzian-shaped voltage pulses to the contact. 

Correlations between electrons in different outgoing 
channels also show up in the distributions of waiting 
times between detections in different channels. We have 
defined multi-channel and cross-channel WTDs and illus¬ 
trated these concepts for a QPC in a chiral setup where 
electrons are injected in either one or both incoming 


channels. In a fermionic Hong-Ou-Mandel experiment, 
where electrons interfere on the QPC, the two-channel 
WTD provides information about the synchronization of 
the sources, while the cross-channel WTD shows signa¬ 
tures of the scatterer and the fermionic antibunching. 

Our formalism relies on the ability to detect single elec¬ 
trons above the Fermi level. A quantum theory of such 
a detector is currently being developed. It would also be 
interesting to make use of WTDs to detect entanglement 
in electronic conductors, for instance by formulating a 
set of Leggett-Garg inequalities for the joint WTDs. 
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Appendix A: Computation of the generalized ITP 

To compute the generalized ITP in Eq. (40), we map 
the outgoing operators onto the incoming ones using the 
Floquet scattering matrix 51,60 

oo Ni 

b a (E)= J2 '52 S F, af i( E ’ E n) & fl( E n), ( A1 ) 

n——oo /3=1 

where D is the frequency of the external driving and we 
have defined the energies 

E n = E + nMd, (A2) 
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for integer n. Since the incoming equilibrium particles 
are described by a Slater determinant, we may evalu¬ 
ate the quantum statistical average in Eq. (40) and we 
thereby obtain a determinant formula 30-32 of the form 

n(£f,;.. -;t s No ,t e No ) = det(I — Q(tf• ■ ■ ;t a No ,t e N J). 

Here I is the identity matrix and Q contains the single¬ 
particle matrix elements of the operator Q a . 

The matrix Q has a block form with elements in the 
(incoming) channel space 

[Q{ti,ti,...,t s No ,t e No )] a p = K a p (A3) 

with the K a « being matrices in energy space. In the 
general case of a time-dependent scatterer with many 
channels described by the Floquet scattering matrix in 
Eq. (Al), these matrices have the elements 


N 0 


[K «p]e,e'=Y, EE S^ a (E m ,E)S F ^(E' n ,E') 

7=1 m=-[E/h£ 2J 
n=— \_E' / hQ\ 

x0 (-E)e(-E)[K(q,t°)] Em _ E ', 

(A4) 

where we have defined 

me,,® ( as ) 

7T £j 


In addition, the floor function is denoted as and we 
have discretized 30,32,4 ' the transport window [Ep,Ep + 
eV] into compartments of width k = eV/Af , where M 
is the number of particles. We always consider the limit 
—> oo. The matrix has a block form with respect to the 
incoming channels, while we sum over the indices of the 
outgoing channels. The result in Eq. (A4) generalizes the 
earlier expression for a single channel 31 to many channels. 


Appendix B: Levitons 


Here we briefly discuss the Floquet scattering theory 
for the periodic creation of levitons.' > 42 ,43,61 Applying a 
series of periodic Lorentzian voltage pulses of unit charge, 


OO 

eV(t)= 

j =—oo 


2hT 

(t-t 0 - jT) 2 + r 2 ’ 


(Bl) 


to one of the reservoirs leads to the creation of one clean 
electron excitation per period — a leviton - - without ac¬ 
companying electron-hole pairs. Here, the period of the 
pulses is denoted as T = 2tt/Q and the parameter to can 
be adjusted to shift the pulses in time. In the adiabatic 
limit, where the variation of the voltage is slow com¬ 
pared to the scattering time, the elements of the Floquet 
scattering matrix can be obtained by Fourier transform¬ 
ing the time dependent scattering phase in Ecp (28) with 
V(t) defined in Eq. (Bl). We then obtain 61 


{ —2e _nflr sinh(fir) if n > 0 
e-° r if n = 0 . (B2) 

0 otherwise 


This Floquet scattering matrix relates particles incident 
on the scatterer to equilibrium particles through Eq. (29). 
Together with the scattering matrix of the scatterer itself, 
it yields the Floquet scattering matrix used for the calcu¬ 
lation of the generalized idle time probability in Eq. ( A4) . 
For the Hong-Ou-Mandel calculations, the time delay be¬ 
tween the sources can be tuned by changing the parame¬ 
ter to in one of the matrix elements in the channel space. 
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